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Abstract. We study the plateau of the magnetization curve at M = M s /3 (M s is the 
saturation magnetization) of the 5=1/2 trimerized XX Z spin chain. By examining 
the level crossing of low-lying excitations obtained from the numerical diagonalization, 
we precisely determine the phase boundary between the plateau state and the no- 
plateau state on the A — t plane, where A denotes the XX Z anisotropy and t the 
magnitude of the trimerization. This quantum phase transition is of the Berezinskii- 
Kosterlitz-Thouless type. 
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1. Introduction 

In recent years the quantized plateau of the magnetization curve of spin chains has been 
attracting much attention. Hida JIJ numerically studied the 5 = 1/2 ferromagnetic- 
ferromagnetic-antiferromagnetic trimerized Heisenberg chain and found the plateau 
of the magnetization curve at M = M s /3 (M s is the saturation magnetization) 
for some parameter region of Jf/Jaf, where Jp and Jaf are the ferromagnetic and 
antiferromagnetic couplings, respectively One of the present authors (K.O.) [[| 
analytically investigated Hida's model to clarify the mechanism for the appearance and 
disappearance of the M = M s /3 plateau. Later related numerical and theoretical [§, 
El H> S B HI are reported in the literature. The magnetization plateaus are also found 
experimentally in 5 = 1 Ni compound [Ni 2 (Medpt)2(/i-ox)(/i-N 3 )]C104 • 0.5H 2 O p 
and in 5 = 1/2 Cu compound NH4CUCI3 flOf . The behaviour of the magnetization 
curve of NH 4 CuCl 3 is quite remarkable, because magnetization plateaus observed at 
M = (3/4)M s and M = (1/4)M S but not at M = and M = (1/2)M S . 

Oshikawa, Yamanaka and Affleck gave the necessary condition for the 
appearance of the magnetization plateau 

n(S — (m)) = integer (1) 

where n is the periodicity of the state, S the magnitude of spins and (m) the average 
magnetization per one spin. Since (|1|) is the necessary condition, it depends on the details 
of the models whether the magnetization plateau exists or not, even if the condition ([]]) 
is satisfied. 

In this paper we study the M = M s /3 plateau of the 5 = 1/2 trimerized XX Z 
spin chain described by 

3j-2,3j-l (A) + ^- 2> 3i-i(A)] + J^, 3j -+i(A)} (2) 

3=1 

where 

h hm {A) = S?S* m + SfSy m + . (3) 

Our model is sketched in figure [I]. 

J J' J' J J' 

-o — o — o — o — o — o 

7 (l+20 J (l-t) J (l~t) / (l+2?) J Q (l~t) 

Figure 1. Sketch of the trimerized XXZ chain. The expression of the lower line 
corresponds to the parametrized form (Jij). 

In §2 we qualitatively discuss the properties of the transition between the plateau 
and no-plateau states by use of the bosonized Hamiltonian. In §3 we determine the 
phase boundary of between the plateau and no-plateau states from the numerical 
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diagonalization data by examining the crossings of the low-lying excitations [|11[|. §4 
is devoted to discussion. 

2. Transition between the plateau state and the no-plateau state 

It is convenient to parametrize the Hamiltonian (||]) as 

L 

H = J {(1 - *)[V- 2 ,3i-l(A) + 2,3j l(A)] 

+(l + 2t)/i 3i)3i+1 (A)}, (4) 

where 

2/ + J f — J 

Jo = ^~ t = -2fTj- (5) 

The model is sketched in figure |I[ The bosonized expression of the Hamiltonian (H) can 
be obtained by the following procedure: 

(a) Transforming (^) into the spinless fermion expression by use of the Jordan- Wigner 
transformation. The spacing between the neighboring spins is taken as the unit 
length. 

(b) Linearizing the dispersion relation of the spinless fermions u(k) = Jq cos A; around 
k = ±&f, where kp = it/ 3 corresponds to the band filling of M — M s /3. The Fermi 
velocity at k = kp is Vp = (\^3/2)J . 

(c) Taking the effects of trimerization and the interactions between fermions into 
account through the procedure similar to that of the standard bosonization technic. 

From the above procedure, we obtain the following sine-Gordon Hamiltonian 



H 

2tt 



where v s is the spin wave velocity of the system, n is the momentum density conjugate 
to 0, [4>(x),H(x')] = i5(x — x'), and the coefficients v B , K, and are related to J , t 
and A as 



where 




VZJv^fAC K = —J- y<f) v s = 27TJ t (7) 



A = — 1 + ^A C = 2tt 1 - ^A . (8) 

The dual field 6 is defined by d x 6 = nil, and we make the identification </> = </> + ^/2n, 
9 = 9 + \/2n. We note that the umklapp term (which exists in M = case and is 
important to describe the transition between the spin-fluid state and the Neel state) 
does not exist, because 2kp is not equal to the reciprocal lattice wave numbers. The 
field (ft is related to the fast varying (in space) part of the spin density S z (x) in the 
comtinuum picture as 



~{cos(Wr-| + V20) +^} (9) 
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which makes it clear the physical meaning of 0. We note that the slowly varying part 
of the spin density is proportional to d<p/dx. 

As is well known, the excitation spectrum of the sine-Gordon model is either 
massive or massless depending on the values of K and y^. In the massive case the M s /3 
magnetization plateau exists, and in the massless case it does not 0. It is convenient 
to discuss the properties of (^j) in the framework of the renormalization group method. 
The renormalization group equations for (P) are 



dK(L)- 1 



dyd L ) 



K{L)' 



y<t>( L ) 



(10) 



dlnL 8* vy ' dlnL \ 2 

where L is an infrared cutoff. Denoting K(L) = 4(1 + y (L)/2) near K(L) = 4, we 
obtain 



dy (L) 



-y<t>i L Y 



&y<t>{L) 



-yo{L)y^{L) 



(11) 



dlnL avy ' dlnL 
and show its flow diagram in figure [| The Berezinskii-Kosterlitz-Thouless (BKT) 




Figure 2. Flow diagram of the renormalization group equation (|ll|). The thick solid 
lines show the BKT lines and the thick dotted line the Gaussian line. 

transition occurs at y§ = \y^\, shown by thick solid lines. At the BKT transition point, 
by substituting y = \y,p\ into eqaution (|TT|), we have 

yo(L) = (T y ° (12) 

y ln(L/L ) + 1 

where y is the bare value. When y < (i.e., K < 4), any small (but not equal 
to zero) amount of trimerization brings about the magnetization plateau. The phase 
boundary between two plateau regions is Gaussian line (thick dotted line), on which 
the critical exponents vary continuously. In the no-plateau region, the effect of the 
trimerization vanishes in the sense of the renormalization group due to the strong 
quantum fluctuations. 

We note that it is dangerous to apply the conventional phenomenological 
renormalization group method to the BKT transition, as is fully discussed in ||12|| . 
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3. Numerical approach 

The scaling dimension of the primary field O m ^ n = exp(m-y/20 + n\[2&) for = is 
given by 

x n ,m = Trm 2 + -^n 2 (13) 
Z ZK 

where n and m are integers with the periodic boundary condition (PBC). According to 
the finite size scaling theory by Cardy ||T3| , [14]], the excitation energy of the finite size 
system at a critical point is related to the scaling dimension as 

Xm,n{L) = ^(E m>n (L) - E g (L)) (14) 

where E g (L) is the ground state energy of L-spin system with PBC. Near the BKT 
transition (K 4), the excitation energy is written as 

-^AE mfi (L) = 2m 2 + y (L)m 2 (15) 
Ztiv s 

— ^— AE 0n (L) = -n 2 - y (L)—n 2 (16) 

2nv s °' V ; 8 y0V J 16 K J 

for integer m, n. Thus, considering equation (|T2|) , we have the logarithmic corrections 
for finite size spectrum. 

To determine the BKT transition point, we use the method developed by Nomura 
and Kitazawa JTl|, in which the level crossings for some excitations are used. With the 
twisted boundary condition (TBC) S^f+i = —Si' y , S^ L+1 = Sf, the integer m in the 
operator C mi „ shifts to m + 1/2 as O m ^ n — > O m +i/2,n- For the scaling dimensions of the 
operators a/2 cos(0/a/2) and a/2 sin(0/y / 2) we have the following finite size corrections 

1 1 , n 1 

2 + 4^ (L) + 2 J 

(17) 

11 1 

^1/2,0 ( L ) = 2 + 4^o( L ) _ ' 

Note that scaling dimensions Xy 20 are not the form ([15]). This comes from the first 
order perturbation of the second term (cos a/20 term) in equation ([]). Denoting 
y,p = ±?/o(l + w) where w measures the distance from the BKT transition point, we 
have for y^ > 

13 / 2 

X \/2A L ) = n + 7Vo{L) 1 + -W 



X l/2fli L ) = n + -Vo( L ) + nV<f>( L ) 



(18) 

1 1 



^i / 2,o( L ) = r4 2/o(L)(1+2w) 



and for y^ < 

1/2,0 = g ~ 4 



xl n „{L) = \--v {L){l + 2w) 



X l/2,0 



(19) 

(£) = ~ + lvo(L) (l + Hu,) . 
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On the other hand, from equation flT6| ) the scaling dimension of £?o,±2 is given by 

xoML) = \~ \yo(L) (20) 

from which we see that £o,±2 and x c {j 2 o ( s f° r V<t> > and c f° r 2/<£ < 0) cross linearly at 
the transition point [w = 0). 

In order to identify the excitation with those of the sine-Gordon model @, we can 
use the following symmetry. The Hamiltonian with PBC is invariant under the spin 
rotation around the S z axis, the translation by three sites, (Sj — > Sj + s), and space 
inversion (Sj — > Si-j+i)- Corresponding eigenvalues are M, the wave number q, and 
P = ±1. The space inversion in the sine-Gordon model are 

0^-0 e^e + 7i/V2 x^-x. (21) 

The magnetization M is related to n as n = M s /3 — M. The "ground state" energy E g 
is the lowest one with [M = M s /3, q = 0,P = l]. 

In our model, the energy level corresponding to the operator Oo,±2 is Eq(M s /3 ± 
2,0,1), where E (M,q, P) is the lowest energy with [M, g, P]. However, we cannot 
directly compare the energies with different M. In the language of the spinless fermions, 
the difference in M corresponds to the difference in the number of fermions, N. Thus to 
compare the energies with different M, we should use E — fiN, where /i is the chemical 
potential of the spinless fermions. Since /i near M s /3 is expressed as 

,= i{ £ „(f + 2,0,l)-*(£-2,0,l)} (22) 
the excitation energy corresponding to Oq <2 is 



A£ ,2 = {^o (jf + 2, 0, l) - E 0, l) - 2/*} 



= ^ {^o + 2, 0, l) + £ - 2, 0, l) } - Eo 0, l) . (23) 

Just the same expression is obtained for C\),-2- Equation (123|) can be also obtained by 
use of the Legendre transformation E — > — if M . 

The excitation energies corresponding to the operators v2 cos(0/v / 2) and 
v2sin(0/\/2) are obtained by the two lowest energies AE(M,P) with the twisted 
boundary condition as 

AE^W^.l) -2^,0,1 



i/ 2 ,o - - ^ 3 > - ^ 3 

(24) 

A^ A0 = ^ TBC (f ,-l)-^(f,0,l 

where E(M S , 0, 1) is the lowest energy with PBC. The excitation energies AEq } ± 2 and 
A-E1/2 o ( s ^ or 2/0 > and c for < 0) should cross linearly at the BKT transition 
point. 

Figure |3| shows the behaviour of AE 0i ± 2 and AE% , 2 for L = 18 spins as functions 
of anisotropy parameter A when t = —0.25. From the crossing point, we obtain 
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Figure 3. Ai?o,±2 and AE^, 2 for L = 18 spins as functions of anisotropy parameter 
A when t = —0.25. From the crossing point we obtain A C (L = 18) = —0.8389. 
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Figure 4. Extrapolation of A c to L = oo when t = —0.25. We obtain A c 



-0.839. 



A c = —0.8389 for L — 18 spins. The BKT transition point for the infinite system 
can be obtained by extrapolating the A c data to L = oo, as shown in figure f|. Thus, 
we can obtain the phase diagram on the A — t plane as shown in figure ^. The point 
M (A = —0.729) is the multicritical point where two BKT line meet together into 
the Gaussian line (shown by thick dotted line) on which the critical exponents vary 
continuously. The A < — 1 region is the ferromagnetic region. 

Let us confirm the conformal anomaly c = 1 which is related to the leading finite 



size correction of the "ground state" energy with PBC as [|Tg, [R; 

7Tf s C 



E S (L) = Le e 



QL 



(25) 
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Figure 5. Phase diagram on the A — t plane. Closed circles are the BKT transition 
points determined from the numerical data as explained in the text. Thick dotted line 
denotes the Gaussian line. The multicritical point M corresponds to the point O of 
figure |2[ 



where e g the energy per one spin for the infinite size system. The spin wave velocity v s 
can be obtained by 

LAE(q = 2ix/L) , N 

v s = lim ^ l—L (26) 

where AE(q = 27r/L) is the lowest excitation energy having the wave number q = 2n/L 
in the M = M s /3 space. Thus we can check the value of c by use of equations (p5j) 
and (0). We have found that c = 1 is realized on the BKT line within the error of a 
few percent. For instance, in case of (t, A) = (0.5, —0.881) on the BKT line, we obtain 
v s c = 0.212 J through equation (]25[) and v s = 0.217 Jo through equation (p6|) . Figure || 
shows the spin wave velocity v s on the BKT transition line. 

From equations (|17[) and (p0|), we can eliminate the leading logarithmic correction 
at the transition points (w = 0) using the following average 

^1/2,0 (-k) + x i/2,o(-^) _ 1 



4 2 

^1/2,0 \L) + Xy 2fi (L) ^ 1 



for y (f) > 

(27) 

for < . 



4 2 

This relation is appropriate to check our method of analyzing the numerical data. The 
averaged scaling dimension equation fl2"7p on the BKT line is shown in figure [JJ We 
can see that the averaged scaling dimension is very close to 1/2, which guarantees the 
consistency of our numerical method. 
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Figure 6. Spin wave velocity on the BKT line. 
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Figure 7. Scaling dimension on the BKT line. Closed circles are x^^ 20 (L), open 
circles x^ 2 (L), and closed squares the averaged scaling dimension 



4. Discussion 

We have obtained the phase diagram on the A — t plane as shown in figure [5]. Two 
BKT lines meet together into the Gaussian line at the multicritical point M where 
A = Am = —0.729. We can analytically predict the value of Am from K = 4 with 
equations (0) and (|8|). We analytically obtain A M = — 3\/37r/21 = —0.777, which shows 
fairly good agreement with the numerical value . The analytically predicted value of the 
spin wave velocity at the multicritical point is obtained by substituting Am = — 3v^3tt/21 
into equation (0), which results in t> s = 0.494J . 

The slopes of the BKT lines for t > and t < near the multicritical point M are 
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the same with each other. This can be explained from the symmetry of the bosonized 
Hamiltonian (||). Hamiltonian (|6|) is invariant under the transformation t <-> —t and 
y/2(j) <-> y/2(j) + n. As the multicritical point M is gone away, on the other hand, the 
BKT lines on the upper and lower planes are asymmetric with each other, as can be seen 
from figure |^. This is quite reasonable because the t <->• —t symmetry does not hold in 
the original spin Hamiltonian (H). From the standpoint of the bosonized Hamiltonian, 
this comes from the existence of higher order terms |17| cos(2y / 20), cos(4v^2</>), ■ • •, of 



which coefficients are also proportional to the trimerization parameter t. If these higher 
order terms are taken into account, the symmetry of the bosonized Hamiltonian under 
the transformation t <-> —t and v2</> <-> y/2(j) + 77 is lost, which explains the asymmetry 
of the BKT lines. 

The mass-generating term cos(v / 2</>) in the bosonized Hamiltonian (||) comes from 
the J x — J x and J y — J y couplings of the trimerization. Strictly speaking, there exists 
another mass-generating term 

— / dx(V0) 2 cos(V20) (28) 

71 J 

which comes from the J z — J z coupling of the trimerization. The effects of equation (p8|) 
and the cos(a/2</>) term in equation (|(J) are mutually competing when A < 0. In 
the 5 = 1/2 ferromagnetic-antiferromagnetic alternating chain ]Tj| [19], this kind of 



competition brings about the transition between the Haldane state and the large-D 
state. In our case, however, the term of equation ( PB[ ) only works to reduce the coefficient 
of cos(v / 20) in equation (j^j. Most simple treatment may be the approximation 
(V0) 2 cos(v^2</>) => ((V</>) 2 ) cos(a/20). In fact, we obtain the phase diagram figure || 
for the "x?/-trimerization model" in which the trimerization exists only in the J x — J x 
and J y — J v couplings and not in the J z — J z coupling. We see that the trimerization 
effect is reduced in figure || in comparison with figure [|, because the no-plateau region is 
wider in figure [| than in figure |8|. We note that the phase boundary of the ferromagnetic 
region is no longer A = — 1, because the SU(2) symmetry is broken even at A = — 1 
in the x?/-trimerization model. This situation is similar to the 5=1/2 XX Z chain 



under the staggered magnetic field p0| . This phase boundary can be calculated from 
the instability of the ferromagnetic state against the M = M s — 1 spin wave excitation, 
resulting in 
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